Soft Ellipsoid Model for Gaussian Polymer Chains 
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A soft ellipsoid model for Gaussian polymer chains is studied, following an idea proposed by Murat 
and Kremer [J. Chem. Phys. 108, 4340 (1998)]. In this model chain molecules are mapped onto 
ellipsoids with certain shapes, and to each shape a monomer density is assigned. In the first part of 
the work, the probabilities for the shapes and the associated monomer densities are studied in detail 
for Gaussian chains. Both quantities are expressed in terms of simple approximate formulae. The 
free energy of a system composed of many ellipsoids is given by an intramolecular part accounting for 
the internal degrees of freedom and an intermolecular part following from pair interactions between 
the monomer densities. Structural and kinetic properties of both homogeneous systems and binary 
mixtures are subsequently studied by Monte-Carlo simulations. It is shown that the model provides 
a powerful phenomenological approach for investigating polymeric systems on semi-macroscopic time 
and length scales. 

PACS numbers: 



I. INTRODUCTION 

For many applications it is important to understand 
polymer dynamics on semi-macroscopic time scales cor- 
responding to configurational changes of polymer chains 
on lengths scales comparable with or larger than the ra- 
dius of gyration Rq- A problem of active current research 
in nano-technology, for example, is the tailoring of thin 
polymer films on surfaces jl], ^ . Spontaneous phase sep- 
aration processes of incompatible polymer blends may be 
used to translate a chemical pattern on the surface (as 
e.g. produced by the micro-contact printing technique) 
into a pattern of varying polymer compositions || . 

Numerical investigation of such problems by means 
of molecular dynamics simulations of semi-microscopic 
bead-spring models is not feasible. An overall configura- 
tional change of an entangled polymer chain composed 
of typically N — 10 3 monomers takes a time of order 
N 4 t v « Is, where t v w lps is a characteristic time for 
vibrational atomic motions. Changes of the morphology 
(e.g. significant changes of phase domains during demix- 
ing of incompatible polymer blends) then take several 
minutes, and this is far beyond accessible time scales in 
computer simulations. 

To overcome this problem, one needs to study coarse- 
grained models. For describing spontaneous phase sep- 
aration processes in polymer blends the Cahn-Hilliard 
equation is often used. In more complex situations in- 
volving hydrodynamic or elastic interactions, generalized 
time-dependent Ginzburg-Landau equations may be em- 
ployed. Exact analytical treatments based on these equa- 
tions are possible for the linear short-time regime and can 
become rather rich for thin films Q . On the other hand, 
the exploration of the large-scale dynamics by numer- 
ical means is hard even within Ginzburg-Landau type 
treatments ||. Nevertheless, numerous theoretical stud- 
ies have been carried out in the past, in particular for 
investigating the behavior in confined geometries (for a 
recent review on this subject, see (6|). 



An alternative way to describe polymer dynamics on 
semi-macroscopic time scales was suggested by Murat 
and Kremer 0. They proposed a model, in which each 
polymer is described by an ellipsoid that can change its 
shape, position and orientation. The probability for a 
particular conformation follows from an intramolecular 
free energy functional. A continuous monomer distribu- 
tion is assigned to each ellipsoid for a given conformation. 
The ellipsoids can penetrate each other and their inter- 
action energy results from a pair potential. If the pair 
potential is given by a delta function, the interaction en- 
ergy between two ellipsoids is proportional to the overlap 
integral of the respective monomer densities times an in- 
teraction parameter e 0. A nice feature of the ellipsoid 
model is that the input parameters (cf. Sec. Ill) can be 
determined from more microscopic models appropriate 
to describe the polymer system on short time scales. In 
fact, in their original work Murat and Kremer consid- 
ered a bead spring model as a starting point and showed 
how the kinetics can be explored on long time scales by 
employing the associated ellipsoid model. 

In the present work we will present a simpler ansatz 
by considering Gaussian chain molecules as underlying 
entities. We will show that the associated "Gaussian el- 
lipsoid model" (GEM) provides a powerful phenomeno- 
logical approach for describing the long-time kinetics of 
interacting macromolecular systems. Moreover, Gaus- 
sian chains serve as a reference system for more real- 
istic microscopic "input models". Compared to these, 
the GEM has the following advantage: By including self- 
interactions of the ellipsoids (see Sec. IV), single ellip- 
soids show Flory type behavior, Rq ~ iV 3 / 5 , while in 
dense systems with monomer concentration c larger than 
the overlap concentration c* they exhibit a random walk 
type behavior, Rq~ N 1 / 2 . In contrast to the bead spring 
model considered in the Rq^N 1 / 2 scaling in dense 
systems is obtained without an appropriate tuning of the 
interaction parameter e with the chain length TV. In 
fact, we will show that the GEM obeys the well-known 
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scaling relation Rg{c, N)~ 7V 3 / 5 /(ciV 4 / 5 ) describing the 
dependence of Rq on the polymerization degree N and 
monomer concentration c || . In this way the GEM allows 
one to study diffusive dynamics on large length scales, 
while accounting for the essential physical properties of 
polymeric systems. 

From a technical point of view, we will show that the 
input quantities for the GEM can be described by explicit 
approximate formulae (valid for iV>30). This allows us 
to calculate the overlap integrals of monomer densities 
analytically, thus speeding up the simulation consider- 
ably. In addition, the explicit form of the input quan- 
tities is helpful for employing the GEM in applications 
and may allow for analytical calculations in the future. 



II. SOFT ELLIPSOID MODEL 

The basic idea of the soft-ellipsoid-model is to map 
one polymer onto one soft particle, that is to do an ad- 
ditional coarse-graining of a microscopic model for poly- 
mers. Soft means, that one ellipsoid can overlap strongly 
with another. With respect to the setup of the model 
outlined in this section, and the definitions and nota- 
tions involved, we follow closely the work of Murat and 
Kremer j7|. 

Given a polymer with (JV + 1) monomers at posi- 
tions y' fe ' in a laboratory fixed coordinate system, k = 
0, . . . ,N, the associated ellipsoid has the center of mass 
r = Ylk=o I (N+l) an d its spatial extension is deter- 
mined by the eigenvalues S = (Si ,82,83) of the radius of 
gyration tensor 



S, 



a(3 



where is the a'th component of y( fc ). Throughout 
this paper we denote spatial indices by Greek letters and 
particle indices by Latin letters. The orientation of the 
associated ellipsoid in space is given by the orientation of 
the principal axis of S a p . Denoting the corresponding ro- 
tation matrix by 1Z, the transformation from coordinates 
y of the laboratory fixed coordinate system to coordi- 
nates x of the ellipsoid's coordinate system is 



1 N 
N + l ^ Kya 



(1) 



fc=0 



x = Tl (y - r) . 
Without loss of generality we choose 
Si > 02 > 5*3 , 



(2) 



(3) 



that means we order the eigenvalues so that the first co- 
ordinate in the ellipsoid's coordinate system always refers 
to the largest axis, while the third coordinate refers to 
the smallest axis. 

Finally, a function g(x; S) is assigned to each ellip- 
soid in its coordinate system, which specifies the mean 



monomer density of the polymer for given eigenvalues S. 
This means, that two ellipsoids i,j with Si = Sj have 
monomer densities g(x;Sj), g(x;Sj) that can be trans- 
formed into each other by a spatial rotation. In the 
laboratory fixed coordinate system the monomer density 
g'iiy) of polymer i is 



0j(y) = e(fti(y-ri);SO 



(4) 



The probability P(S) for a polymer to have eigenvalues 
S and the conditional monomer density p(x; S) are the in- 
put quantities for the soft-ellipsoid-model. They have to 
be determined from a microscopic model and their form 
for the GEM (Gaussian chains) is discussed in Sec. [II. 

The free energy functional of an ensemble of M soft 
particles in a volume V is divided into an intramolecular 
part, which accounts for the possible internal configura- 
tions, and an intermolecular part, which describes the 
interaction between the polymers, 

F = -Fintra + -Winter ■ (5) 

The intramolecular part is given by the probability P(S) 



A I 



A I 



^intra = £ i^La = ~k B T £ In P(S() , (6) 

i=l 1=1 

where ks is the Boltzmann constant and T is the tem- 
perature. The intermolecular part is given by 



1 M •■) (l M 

Pinter = 2 X! Writer + I o ^ 

i=l (j^i) \ i=l 



F 



(a) 



F 



0¥0 



d A y d 6 z v(y,z) e ' l (y)g' J {z) 



(8) 



We wrote "(J 7^ £)" in eq. ([?]), since one may or may 
not include self-interaction terms. The effective pair po- 
tential v(y, z) between two monomers can include indi- 
rect contributions mediated by solvent molecules. For 
the most simple choice accounting for excluded volume 
effects ||, v(y,z) = eb 3 S(y — z) with e and 6 3 being a 
"contact energy" and a "contact volume", respectively, 
the interaction is given by the overlap integral of the 
monomer densities, 



F 



(y) 



'y ei(y)ej(y) 



(9) 



Equations fix the thermodynamics of the model. 

Generalizations to polymer mixtures are straightfor- 
ward 0. In a binary polymer blend we have two types 
of polymers A and B. For simplicity we will assume 
that both types have the same polymerization degree, 
Na = Nb = N, and that a polymer of each type inter- 
acts with a polymer of the same type with an interaction 
tAA = £BB = t, while polymers of different types interact 
with cab = e{l + 6) (see Sec. |v|). 

To summarize, the external parameters defining the 
thermodynamic state of a homogenous system are the 
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FIG. 1: Numerical results of the scaled probability function 
Pr(Rq) for N = 30(D), 100(B), 300(o) and 1000(»). The 

approximate function Pr(Rg) defined in the text (cf. eq. ( |l3| ) 
and Tab. |) is drawn as a straight line. The inset shows the 
same data in semi-logarithmic form to prove the validity in 
the range of very small and very large R G . 



polymerization degree N, the overall monomer concen- 
tration c = (N + 1)M/V, and the reduced interaction 
strength e = e/ksT. In a binary blend we in addition 
have the mismatch interaction S and the fraction j a (or 
concentration ca) of A polymers. 

For describing kinetic properties we implement a dis- 
crete time Monte-Carlo algorithm involving three dif- 
ferent types of moves of a randomly chosen ellipsoid 
j/J: (i) A translation, where the center of mass is dis- 
placed by a random vector Ar, the components of which 
are drawn from a uniform distribution in the inter- 
val [— Ar max /2, Ar max /2]; (ii) A rotation of the ellip- 
soid with equal probability in the steradian 47r; (Hi) 
A change of the ellipsoid's size, where a triplet of ran- 
dom numbers AS a uniformly distributed in the inter- 
val [max(— A5' max . ct /2, — S a ), S' maXiQ /2] is added to the 
eigenvalues S. We restricted the changes AS a to values 
that result only in positive S a . Moreover, for a move 
according to rule (Hi) it is possible that the order of the 
principal axes changes because of condition (^), which 
requires a simultaneous update of the orientation of the 
ellipsoid. All possible moves are attempted with the same 
probability and we define one Monte-Carlo step (MCS) 
as 3 M trials to change the state of an ellipsoid. 

We note that the dynamics are purely diffusive and 
neglect advective processes due to hydrodynamic flows. 
Moreover, the dynamics do not capture entanglement ef- 
fects, which are important on time scales smaller than the 
disengagement time td needed for a polymer to diffuse 
over a distance of order Rq M . As a consequence, when 
adjusting the Monte-Carlo time step to real experimen- 
tal time scales, one should take care that the dependence 
of td on N will generally be different in the simulations 
and the real system. 



III. GAUSSIAN ELLIPSOID MODEL (GEM) 

The Gaussian chain is one of the best examined models 
for polymer chains. In this model the distribution of bond 
vectors Ay between Kuhn segments is given by 

^=vS exp (~^)' < 10 > 

where for simplicity we have identified the root mean 
square bond length b with the parameter specifying the 
contact volume in eq. (||). We define b as our length unit, 
i.e. 6=1. 

In the following we discuss the input quantities 
P(S,N) and g(x;S,N). For clarity we have now ex- 
plicitly marked the dependence on N, but we take the 
freedom to suppress the argument N, whereever it seems 
to be appropriate. Explicit approximate formulae are de- 
termined for the input quantities based on Monte Carlo 
simulations, as well as asymptotic expansions and numer- 
ical results published in earlier work (see below). The 
simulations are conducted following the procedure de- 
scribed in [ pTj[ . Averages were performed over typically 
10 7 different chain conformations. 



A. Probability distribution P(S) 

The exact analytical form of the distribution function 
P(S, N) for Gaussian chains is rather complicated. It 
can be presented in terms of multiple integrals , but 
an exact calculation from the respective formulae is very 
cumbersome from the numerical point of view. Useful 
explicit expressions are available only in the limit of large 
eigenvalues S a jfl], [l2| . To obtain a convenient formula 
for practical computations, we first consider the simpler 
distribution function Pr(R g , N) of the squared radius of 
gyration 

R 2 G = tr S af} = S 1 + S 2 + S 3 , (11) 

which can be expressed in terms of a single integral [fl3"[ . 
For large N — ► oo, Pr(R g , N) obeys the random walk 
scaling 

/>oo 

P R (R 2 G ,N)~N- 1 P R (B$;/N), / duP R {u) = l. (12) 

Jo 

The scaling function Pr(u) behaves as Pr(u) ~ 
u"° exp(— Co/u) for u — > and Pr(u) ~ u v °° exp(— u) 

TABLE I: Parameter defining Pr(Rg) in eq. @ 



tir 


1.0 


dR 


0.08020 


d R 


1.842 


K (2d R ) 


0.015923 
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for u^co O, 15 with Co, Coo, vq, and being con- 
stants. Note that the first moment of P R {u) must equal 
1/6 to obtain the well known result Rq = {Rq) 1 ^ 2 ~ 
(7V/6) 1/2 in the limit of large N. 

In this section, (...) denotes an average over an ensem- 
ble of Gaussian chains, while in Sees. |v|,|v| it will also be 
used with respect to an ensemble of interacting ellipsoids. 
From the context it will be clear, what kind of average is 
meant. 

Since the exponential terms dominate the behavior of 
Pr(u) for small and large u we make the ansatz 



a) 



Pr(u) = 



riR-l 



u-"* (a R d R ) 
2K rlR _ 1 (2d R ) 



■ exp 



u 2 a R 

"r — 

a R u 



(13) 



where K nn _i{.) denotes the modified Bessel function of 
(n R — l)th order accounting for proper normalization. For 
simplicity we restrict n R to half-integers, but it is not 
supposed to equal — vq or —Poo- Since the form of the 
maximum turns out to depend sensitively on n R , we de- 
termined n R by a least square fit of (Oft to the Monte 
Carlo data in the region given by the full width at half 
maximum. After determining n R = 1 by this procedure, 
we calculated the constants a R and d R from the first two 
moments of P R (u) (see Appendix A and Tab. |jj). 

Fig. [I] shows that the scaling (|I^) is well obeyed for 
N > 30 JT| and that eq. |[|) provides an excellent ap- 
proximation to the simulated data. Moreover, higher mo- 
ments (Rfj 1 ), m = 2, 3, . . . calculated by using eq. Jl3| ) 
deviate by less than one per cent from the exact result 
(see Appendix A). 

To find a useful approximate formula for P(S, N), we 
first consider the distributions 



P a (S a ,N) = J d 2 S p P(S,N) , 



(14) 



of the single eigenvalues S a , a = 1,2, 3. In view of eq. < \l l| ) 
these can be expected to behave similar to P^R^jN). 
Using the random walk scaling for large N, we write 
P a (S a ,N)~N- l P a (S a /N), f™duP a (u) = l, with P a (.) 
as in eq. (jig), except for the fact that the parameters n R , 
a R , and d R are replaced by n a , a Q , and d a . These param- 
eters have been determined analogous to the procedure 
described above (for details see Appendix A) and they 
are listed in Tab. [n| Again, there is an excellent agree- 
ment between these functions and the simulated data, 
see Fig.@. 
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FIG. 2: Comparison of Monte-Carlo results for the scaled 
probability functions P a (S a ) for N = 30(D), 100(B), 300(o) 
and 1000(») with the approximate formulae defined in eq. (113) 
and Tab. 0. The insets show the respective data in semi- 
logarithmic form again. 



TABLE II: Parameter defining P a (S a ) in eq. (|l|) 



Sa 


Si 


s 2 


S 3 


n a 


1/2 


5/2 


4 


a a 


0.08065 


0.01813 


0.006031 


d a 


1.096 


1.998 


2.684 


Jf„ a _ 1 (2da) 


0.094551 


0.0144146 


0.0052767 



Next, in order to obtain the combined probability den- 
sity P(S, N), we neglect correlations among the S a , and 
write 

p(s, n) = n p a (s a , n) = n i p a (is) 
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FIG. 3: Scaled contour plots showing the comparison of 
\n(P a (Sa)Pp(Sp)N 2 ) (straight lines) with numerical data of 
In P(S a , 5/3 ) (dashed lines). The isolines are drawn at func- 
tion values of 10 n , where n is an arbitrary integer number. 
The agreement of the separation ansatz eq. ( |l5|) with the nu- 
merical data is quite good especially in the region of large val- 
ues, which have the highest statistical weight. However the 
separation ansatz neglects condition eq.(^|). For large values 
of Sa the Monte-Carlo data become noisy due to insufficient 
statistics. 



A useful test of the accuracy of eq. (|15| ) is to check 
whether the correlation coefficients 



(5) _ (SgSp) - (S a )(Sp) 

- (S a S ) 



(16) 



are much smaller than one. The results summarized in 
Tab. Ill show that all are less than 6%. We also 



checked whether the asphericity |l7j 
(R 2 G ) - 3(S 1 S 2 + 



(A) = 



S?S; 



2^3/ 



(R 2 



(17) 



Gl 



of the ellipsoids is significantly changed compared to the 
true value. From the simulations we find (A) = 0.527, 
while from eq. ( |l5|) it follows (A) =0.542, which amounts 
to a relative error of about 3%. 

A more sensitive test is the direct comparison of the 
combined probability density 



P(S a ,S ) 



dSy P(S) 



77^ a, j3 



(18) 



with P a {S a )Pi3(Sp). As shown in Fig. ||, the agreement 
between the ansatz ( |l5| ) and the true behavior obtained 
from the simulations remains good, except for the regions 
where P(S a ,Sp) becomes extremely small. We plotted 
the logarithm of P(S Q ,S» in Fig. |, since we are in- 
terested in determining Fi n tra cx lnP(S). The main re- 
striction of ( |l5"| ) is that it does not obey the condition 
@, i.e. P(S) = for S 2 ,S 3 < X and S 3 < S 2 . However, 
this poses no severe problem, because the difference in 
the three scales^ (Si) : (S 2 ) ■ (S3) = 12.07 : 2.717 : 1, 
cf. 



eg. Tab. [TV, 
imum of P(S) 
eq. (jl^) breaks down. 
P a cx exp(— C a /S a ) for small S a , P is very small in the 
region, where it should be exactly zero. 



ensures that the location of the max- 
is well separated from the region, where 
Due to the exponential decrease 



B. Monomer density g(x; S) 

For p(x;S,7V) no exact analytical expression is avail- 
able. Since £>(x; S, N) is a function of many variables, 
various simplifications are necessary to allow for an ef- 
ficient computation of the interaction terms in eq. (|^). 
Murat and Kremer [M proposed a self-similar monomer 

distribution (R a = Sa ), 



e(x ; s,JV) = 1 ^y-e(x/R) ) 



(19) 



where we introduced the short-hand notation g(x/R) = 
q(x\/R\, x 2 /R 2 , X3/R3). In contrast to Q, we included 



TABLE III: Results for the correlation coefficients to test 



)■ 


«SlS2) 

((S 2 S 3 ) 

((8 3 Si) 


-{S 1 ){S 2 ))/{S 1 S 2 ) 
-(S2}{S 3 ))/{S 2 S 3 ) 
-(S 3 )(S 1 ))/(S 3 S 1 ) 


0.054 
0.059 
0.017 
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the factor [(AT+1) / Y[ a Ra] in eq. ( |19| ) due to the normal- 
ization 



J d 3 xg{x;S,N) = N + l 



(20) 



While the self-similarity assumption seems to be 
plausible, we will next show that it is not strictly 
valid. To this end we consider the function 
lY[ a R a ] £>(x; S, N)/(N + 1) in terms of the scaled vari- 
able u = x/R, i.e. we consider 



g(u;S,N)i 



URg 

N + l 



g(R lUl ,R 2 u 2 ,R 3 u 3 ;S,N). (21) 



If the scaling (|n]) would hold true, this function would 
equal /5(u), and hence any averaging over regions in S- 
space would leave it invariant. In particular we define 

Qa/3(u a ;S/3,N) = 

/>CO 

IT diifi / TT dSu P(S, N) £(u; S, N)@2) 

'.._t_ Jo ..-La 



J — i 



that means the scaled monomer density with respect to 
one principal axis a under the condition that one ar- 
bitrary eigenvalue corresponding to the axis (3 is given. 
These are nine different functions, and the Monte-Carlo 
results for some of them are shown in Fig. || for N — 1000. 
For other N > 30 we find almost exactly the same func- 
tions, that means the scaling with N is well obeyed. In 
order to obtain the functions displayed in the figure, we 
divided the S'a-axis in eight different intervals, delimited 



at positions Sjjp 



.S 



(9) 



, s 



(!>) 



oo, each of 



them having the same statistical weight 1/8, i.e. 



dSf) P (S(3) 



1.. 



(23) 



and averaged g a p(u a ', Sp) in each of these intervals. 

Figure ^a shows the eight resulting functions for 
£11(1*1; Si, N). By contrast to the self-similarity assump- 
tion, the curves differ significantly from each other. They 



TABLE IV: Moments of the eigenvalues S as determined by 
Monte-Carlo simulation. The table displays results for chains 
of length iV = 1000 averaged over 10 7 different realizations. 
The data for the radius of gyration Rq are given for the sake of 
comparison. The deviations from the exact values for N — > 00 
are due to statistical errors and finite chain lengths. 



{Rq}/N 0.16685 


(Rg)/N* 0.035302 


{Si)/N 0.12758 
{S 2 )/N 0.028708 
(S 3 )/N 0.010568 


{S'i)/N* 0.023009 
(SaViV 2 0.0010300 
{Si)/N 2 0.00013143 


(SiS 2 )/N' J 0.0038714 
(S 2 S 3 }/N 2 0.00032239 
{SsS^/N 2 0.0013717 


(S'i)/N 3 0.0055805 
(Si)/N 3 4.6206 10~ 5 

(sD/iv 3 1.9338 nr 6 
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FIG. 4: One-dimensional conditional monomer density 
Qap{x a / R a ; Sp), scaled with Ra ~ VSa, to demonstrate the 
breakdown of condition ( p^| ) in a strict sense. The curves were 
sampled for chain-lengths N = 1000. 



exhibit an unimodal shape only for the smallest Si , while 
for the most likely Si, Qn(ui; Si, N) has a bimodal shape 
that becomes more pronounced for larger Si. A similar 
bimodal shape for Gaussian coils has been found for the 
average monomer density g(x, N) = J d 3 S P(S)g(x; S, AT) 
by Janszen et. al. jl(J], and earlier for a model repre- 
senting polypropylene by Theodorou and Suter jy| . For 
£22(1*2! S2) shown in Fig. |^b, the deviations from scal- 
ing are smaller and are appreciable only for the largest 
S2, while for £33(113; S3) in Fig. ||c the deviations are 
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FIG. 5: Scaled one-dimensional density functions for different 
values of N. The s traight lines denote the approximate forms 
given in eq. (HIBi-c). 



insignificant. We also examined the six other functions 
and found that the deviations from scaling are small in 
all these cases. 

In fact, knowing that the shape of p(x, N) in the in- 
direction is of bimodal type, while in the ^-direction 
it is of unimodal type |10), the breakdown of the self- 
similarity concept is not that unexpected. If scaling 
would hold true, the monomer densities with respect to 
some fixed S should correspond to the average monomer 
density g(x 7 N). However, if one considers the case of 
almost equal values of Si and S2, i.e. Si <L S2, the 



monomer densities along the x±- and a^-direction should 
be the same due to symmetry in contrast to the scaling 
prediction. 

Although the self-similarity assumption is not truely 
valid, it can be used as a good working tool for the follow- 
ing reasons. The monomer densities enter the calculation 
of -Pinter via integrals, cf. eq. For these integrals to 
be correctly calculated, it is important that the size of 
the region, where g(x, S, N) is non-neglible, scales prop- 
erly. From Fig. ^ one can see that this is indeed the case. 
Therefore, for our purposes it is sufficient to approximate 
the various functions shown in Fig. |4^,-c for different Sp 
by their average 



dS/3 Pp{Sp) g a p{u a \ Sp, N) 



(24) 



We marked these averaged monomer densities by a tilde, 
since one can view them as "effective scaling functions" 
representing g(u) from eq. (|l9| ) (after integrating out the 
two w^-coordinates orthogonal to u a ). Note that we sup- 
pressed the iV-argument in the definition (|2~i|), since, as 
already mentioned, the dependence on N predicted by 
eq. ([19]) is valid. This validity can also be inferred di- 
rectly from Fig. §, where the data for N = 30, 100, 300, 
and 1000 can not be distinguished. 

The scaled monomer densities g a (u a ) with u a — x a /R a 
shown in Fig. can be conveniently approximated by a 
superposition of Gaussians 
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(25) 
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The parameters er Q , c±, and u\ were determined by least 
squares fits and are given in Tab. 0. While £i(«x) and 
£3(^3) are quite well approximated by eqs. (25|) and (p7|), 
respectively, eq. (^) does not provide such a good de- 
scription. A better account for 02(^2) can be achieved by 
adding shifted Gaussian functions as in ( p5| ) , and is given 
explicitly in appendix |b|. However, we decided to deal 
with eq. (26), because the better description increases the 
computation time of .Pinter by a factor of four, but does 



TABLE V: N umeri cal results for the constants defining 
g a {u a ) in eqs. ( frlIB a-i 
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not much improve the overall accuracy ofp(x; S). We 
note that the functions g a (u a ) in eqs. fllIIB| a-c) should 
not be confused with the function p(x, N) (or its "com- 
ponents" after integrating out two coordinates) that has 
been considered earlier by Janszen et al. |1Q] . 

In order to specify the multivariate scaled density g(u) 
we employ, as for P(S) before, a separation ansatz. This 
yields for p(x; S, N) from eq. (|l9| ) the final result 



(x;S,JV) = (JV+l)JJ 



g(x a /R a ) 



(28) 



where g(x a /R a ) are taken from eqs. ( |[IIB| a-c). The sepa- 
ration ansatz may be used also without invoking the self- 
similarity assumption in order to reduce the complexity 
of the multivariate density g(x; S, TV). 

To test the validity of the separation ansatz we define 



for each Gaussian chain by v a — (N+l) x ^ 



N 

k=0 



\/R 



the (scaled) mean moduli of the monomer coordinates in 
the principal axis system, and check whether the corre- 
lation coefficients 



(v) _ (vgVp) 

lafl — 



M (vp) 



V a Vp) 



(29) 



are much smaller than one. Note that, because of the 
symmetry g a (u a ) = g a {— u a ), we considered the mo- 
ments of the moduli in 
find | 7q/3 ' 

when comparing the two-variate densities 



As shown in Tab. VI 



' 1 to be less than 5% in the simulations. 



we 
Also, 



g a f}{u a ,Uf}) 



(30) 



with the two-fold products g a (u a )g/3(up) in Fig. [| a sat- 
isfactory agreement is obtained. 

As an ultimate test of the separation ansatz we sam- 
pled the full three- variate scaled monomer density g(u) in 
the simulations. To this end, we subdivided the u = x/R 
space into small boxes of length Au a = 0.07 and deter- 
mined the mean number of monomers in each box by av- 
eraging over Gaussian chains irrespective of their S. The 
relative error between the separation ansatz JT g a (u a ) 
with g a (u a ) from eqs. (lIIBa-c) and the simulated g(u) 
is plotted against g(u) in Fig. [?] for all u. We see that 
the relative error is confined to a narrow band of maxi- 
mal 25% discrepancy except for very small g(u). These 
small g(u), however, are not important in the calculation 
of the overlap integrals in eq. (^). In principle, it would 
be desirable to do the same analysis also with P(S) but 



TABLE VI: Test of the ansatz eq.(gg 
ficients denned in eq. (J29|). 



with correlation coef- 



{(V1V2) - («i)(«2»/(uiV2> -0.048 
({v 2 vs} - <U 2 ) W)/<«2«3) -0.027 
({v 3 vi) - (v 3 )(vi))/(v3Vi) -0-018 
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FIG. 6: Contour plots showing the comparison of 
Qa(x a / Ra)Qp{xf3 / Rp) (straight lines) with numerical data of 
Qapixa/Ra, oo g/ Rfj) (dashed lines). The isolines are drawn at 
function values, which are integer multiples of 0.02. 



since one Gaussian chain generates only one realization 
of S, it is very hard to obtain reasonable statistics in this 
case. 

The monomer density g(x, S, N) specified in eq. J28j) 
allows us to calculate the overlap integrals and thus -Fmter 



from (I 
pendix 



) analytically. This is explained in detail in Ap- 



9 




0.01 0.02 0.03 0.04 0.05 
P(u) 



FIG. 7: The plot shows t he re lative error of the approximate 
function defined in eqs. when compared to 3D nu- 

merical data of g(u) for chains of length TV = 1000, versus 
the function value of g(u). For large function values, that 
give the main contribution to -Pinter, the relative error is quite 
small. The lines denote a 25% error interval. 



C. Self-interactions 

Gaussian chains are by definition non-self-interacting 
in contrast to real polymer chains. In order to obtain 
a physical model, we therefore include a self-interaction 
term in eq. (Q), which amounts to an intermolecular part 
of the form 



1 M 



inter 



l 



M 

2 / j into 



(31) 



with F^Hr given in eq. (0). The factors (1/2) in front 
of both the double and simple sum in (Jfy arise to avoid 
double counting of monomer interactions. 

In principle, it would be more appropriate to include 
the self-interaction term into the intramolecular part in 
eq. (^). In fact, when referring to some kind of real- 
istic model for single polymer chains, as e.g. the bead- 
spring model considered by Murat and Kremer H this 
self-interaction part naturally is entailed in Pjntra- I n our 
case, however, interactions between monomers belonging 
to the same or different chains are treated on the same 
footing and we thus decided to include the self-interaction 
part into eq. (|). The GEM thus specified can be re- 
garded as a Flory type model for polymer systems. 



D. Remarks on the algorithm 

In this section we shortly describe some details of the 
algorithm used throughout the rest of this work. As the 
Monte-Carlo method is well established [[l9) , we just com- 
ment on some specialties in this case. 

The main challenge for the setup of an optimal algo- 
rithm is to calculate the interaction term -Pinter efficiently, 



even when the calculation can be done analytically. As 
one polymer in a melt overlaps with many other ones, 
it is crucial to calculate only those interaction terms, 
which contribute significantly to .Pinter' Therefore, the 
monomer density has to be cutoff at a certain length and 
we choose this cutoff to be x atCU t/R a — 2.2, cf. Fig. ||. 
The neighbors of a polymer are defined as those poly- 
mers, which interpenetrate each other within their cutoff 
shape. Because the ellipsoids are strongly aspherical, a 
simple test of the distance between the centers of mass 
is not efficient enough. 

In order to determine whether two ellipsoids are neigh- 
bors, we introduced a grid of cells, with a spatial dis- 
cretization length given approximately by the smallest 
significant scale, i.e. (R3). To every cell in the grid is as- 
signed a list of polymers, that overlap with the cell, and 
conversely, to every polymer is assigned a list of cells, 
which it covers. Determining the neighbors j of one poly- 
mer i then is simple: One goes through the list of cells 
assigned to polymer i, determines all polymers j that 
are associated with each member of this list and finally 
sorts out multiple occurrences of j. For each polymer we 
keep a list entailing the neighbors with the respective in- 
teraction energies, and the total value of the interaction 
energy with the neighbors. In this way the updating of 
the lists after each Monte-Carlo move becomes simple. 
Since the calculation of Pinter consumes the main part of 
the CPU time, we parallelized this part, thus speeding 
up the simulation by a factor of 1.5 on a two-processor 
machine. 

With the method described here, a fast calculation of 
Pinter is possible. However it should be noted, that the ef- 
ficiency depends on the possibility to integrate the prod- 
uct of two densities analytically. The time for calculation 
depends essentially both on the overall monomer concen- 
tration c and on N (or on Rq), because these quantities 
determine the average number of neighbors of one parti- 
cle. 

For performing one elementary step of the Monte-Carlo 
simulation, where we try to change either the orientation, 
the center of mass position or the eigenvalues S of one 
randomly chosen ellipsoid as explained in Sec. II, we still 
have to specify Ar max and AS max . We choose 



K r R G , AS n 



K s (S a ) 



(32) 



with PT r = 0.25, K s = 0.5 and the (S a ) given in Tab. [TV. 
A trial for a change from an initial state i to a final state 
f is accepted with a probability pi / given by the Glauber 
rule, i.e. p lf = {l/2){l~ta,nh[(F f -F l )/2k B T}), where F t 
and Ff are the free energies of the initial and final state, 
respectively. After each trial the time t is incremented by 
To/3M, where To is a characteristic time scale (1 MCS) 
that has to be adapted to the particular physical situ- 
ation under consideration (see the discussion in Sec. IV 
below) . 

Unless noted otherwise, the initial conditions are as 
follows: The center of mass positions of the ellipsoids are 
randomly distributed within the simulation box. The 
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orientations of the ellipsoids are random, and the eigen- 
values Si are drawn from the distribution P(S, N) for the 
non-interacting Gaussian chains given in eg . ( |l4| ) . The S a 
are sorted afterwards according to rule (|3|). Equilibrium 
quantities are determined after proper thermalization. 

To summ arize this section, eqs. (|l^) and ( p8| ) together 
with ( III B ) define the input quantities for the GEM, and 
one has to tak e care of including self-interactions as ex- 
plained in Sec. Ill C| . The explicit form of the respective 
formulae allows for an easy reproducibility of the results 
to be discussed in the following sections. Despite ap- 
proximate, they capture the essential features found for 
Gaussian chains: The normal random walk scaling with 
N, the strong asphericity of the chains, the bimodal form 
of the monomer density and its rapid decay beyond cut- 
off lengths that scale linearly with R a . We note in pass- 
ing that we performed also simulations for self-avoiding 
chains and found similar features except, of course, for 
the scaling with N. While it is possible to provide more 
accurate formulae for the input quantities by introducing 
more parameters, we are convinced that this is of minor 
importance for the study of structural properties and ki- 
netic phenomena on semi-macroscopic length scales. 
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IV. HOMOGENEOUS SYSTEMS 

The aim of the following section is to show that the 
GEM reproduces the most important semi-macroscopic 
properties of dilute and dense polymeric systems: The 
scaling properties of the radius of gyration R G upon N 
and c, the scaling behavior of the distribution P R (R G , N) 
in a dense homogenous system of interacting ellipsoids, 
and the existence and scaling of the correlation hole. We 
partly follow the analysis proposed by Murat and Kremer 

0- 

As mentioned above, the free energy as given in eqs. 
is some kind of Flory-like description of polymers. 
Accordingly, the gyration radius of a free ellipsoid, that 
only interacts with itself but not with other ellipsoids, 
scales as R 2 G ~ N 6 / 5 §, see Fig. |a. 

In order to obtain the dependence on e also, let us 
redo some Flory type calculation, where for simplicity we 
consider soft spheres instead of soft ellipsoids. Within 
this approximation the monomer density from eq. (|2~F 
would become 



Qr{r:R G ) = 



(33) 



with Q r (u r ) = exp(-u2/2o-2)/(27rcr;?) 3 / 2 , a 2 r = 1/3. Such 
a simplified form of the monomer density has been 
discussed already in Substituting P R (R G ) from 

eqs. (|l2| , |l3| ) into the free energy functional F^(R G ) = 



FIG. 8: a) Mean squared gyration radius of free (c — > 0), but 
self-interacting ellipsoids as a function of N, displayed for 
various values of the interaction strength e. The Flory result 
(R%) oc N 6/5 is shown for comparison. The inset shows the 
dependence on e. The line indicates a power law 6(e) oc e 2//5 . 
b) Crossover scaling from Gauss to Flory type behavior. 



-fc B TlnP fl (i?^) + 47re j^°drr 2 gl(r;R G ) we obtain 
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The most probable value R 2 G results from minimizing 
with respect to R G , 

3/2 
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This equation has the asymptotic solution 
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(35) 



(36) 



m reason- 



Since R G ~ flg , we expect R G 
able agreement with the numerical results for free, self- 
interacting ellipsoids, see Fig. Ha. 
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For small e, one should recover the Gaussian behavior, 
i.e. we expect Rq^N for e<Ce*(-/V). On the other hand, 
the scaling ( |36| ) must hold true for e^>e*(./V). Thus, from 
continuity, we find e*(A) ~ A -1 / 2 . The corresponding 
scaling form Rq(N , e) — N f (e 2 N) with /(-u)^ const, for 
u — > and /(li)^it 1 / 5 for u — > oo is verified in Fig. |^b. 

We now discuss systems composed of many interact- 
ing ellipsoids. As is well known from self-consistent field 
arguments Q, the Flory type scaling ( |36| ) changes to a 
normal random walk scaling if the monomer concentra- 
tion c exceeds the overlap concentration c*. That this is 
indeed the case is shown in Fig. ^. Approaching c* from 
below we can write (cjA)" 1 / 3 ~ i? G (A,e) ~ e 1 / 5 ^/ 5 
(for e>A~ 1 / 2 ), i.e. c* ~ e~ 3 / 5 A~ 4 / 5 , so that more gen- 
erally we expect the scaling form Q 

R 2 G (N,e,c) = e 2 / 5 A 6 / 5 /(ce 3/5 A 4 / 5 ) (37) 

with /(A) -> const, for A<1 and /(A)->A- 1 / 4 for A>1 
to describe the crossover from dilute to dense systems. 
Accordingly, R 2 G ~ e 1 / 4 C - 1 / 4 A for ce 3 / 5 7V 4 / 5 > 1. To 
demonstrate the validity of the overall scaling predicted 
by eq. (|^), we show in Fig. |To| as an example Rq/N 6 / 5 
as a function of cA 4 / 5 for a fixed value e = 1.0. As can 
be seen from the figure, the data for different c and A all 
collapse onto a common master curve. We note that the 
scaling (^tJ) is followed only for e> A" 1 / 2 . For e< A^ 1 / 2 
by contrast, the ellipsoids would always exhibit Gauss 
type behavior. 

We further show in Fig. [Il] that the distribution func- 
tion Pb^Rq, A) for dense systems (c ^ c*) obeys the 
same type of scaling as for non-interacting Gaussian 
chains, eq. P R (R 2 G ,N)^N- 1 P R (R 2 G /N). In partic- 
ular, the ansatz ( filf ) can be used equally well for Pb(u) 
up to a change of the parameters a^, and d R (see Tab. ffl 
and the values given in the caption of Fig. [n]) . 

To complete our study of the static properties of dense 
homogeneous systems, we discuss the properties of the 
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FIG. 9: Squared gyration radius Rq as a function of N in 
a dense (c = 0.85) system of M = 1000 -> 4000 ellipsoids 
(double- logarithmic representation). For various values of 
constant e a linear behavior is observed. 
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FIG. 10: Test of the scaling relation between the squared ra- 
dius of gyration and the monomer concentration for various 
values of N and constant e = 1.0 (double-logarithmic rep- 
resentation). For small concentrations the chains are almost 
free and Rq saturates, corresponding to Flory type behav- 
ior. In the region of high concentration the expected scaling 
relation is almost fulfilled. 
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FIG. 11: Scaling behavior of the distribution function P(R G ) 
in a melt system for N = 30 (*), 50 (■), 100 (•), 200 (A) and 

400 (T). The straight line refers to a function of form eq. (p^[), 
but with changed parameters an — 0.0665 and dn = 3.52. 
The inset shows the same data in semi-logarithmic represen- 
tation. The external parameters are c = 0.85 and e = 1.0. 

so-called correlation hole ||. Therefore we determine 
the mean monomer density of all ellipsoids except the 
i'th one as a function of the distance from the center of 
mass of ellipsoid i, 

Q' dist (r,N) = (j2 Q j(r-ri,N)). (38) 

After performing the average in (|38l ) , this mean monomer 
density of "distinct ellipsoids" is equal for all i and de- 
pends on r — |r| only. For large r > Rq , £?dist( r )^ c ' while 
for small r <C Rg > £?dist ( r ) mus t be smaller than c due to 
the fact that ellipsoid i has been excluded from the sum 
ineq. ©. Since An dr r 2 [g' dist (r, N) - c] = -(A+l), 
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we expect a scaling [g' dist (r, N)—c] ~ N~ 1 / 2 f{rN~ l l 2 ) for 
the correlation hole, in reasonable agreement with the 
simulated results shown in Fig. |l2| . 

We conclude this section with some remarks on the 
dynamical behavior of the system. As can be seen from 
Fig. the time-dependent mean square displacement 
([r»(t) — i"i(0)] 2 ) of an ellipsoid exhibits normal diffusive 
behavior for times t > td , where td is the disengagement 
time, i.e. the average time an ellipsoid needs to diffuse 
over a distance Rq- The short time regime t < td is 
not of interest here, since we do not intend to capture 
the complicated dynamics of polymer systems on these 
time scales. In particular, the diffusion coefficient D = 
limt_ >00 ([rj(t) — Yi({))] 2 ) /Qt will usually not exhibit the 
desired scaling with N, e.g. D^N^ 1 for Rouse chains or 
D^N~ 2 for entangled chains. 



In fact, because D depends on At 



and Ar-n 



is chosen to be some fraction K r of R G ~ N 1 / 2 (see 



Sec. Ill D ), the variation of D with N is influenced by 
the simulation procedure itself, i.e. by the value of K r . 
For the choice K r = 0.25 used in our simulations, D in- 
creases with N (see Fig. [l3|). Conceptually, this is not 
a crucial problem. Since the element ary tim e scale tq in 
the Monte Carlo procedure (see Sec. HID) is arbitrary, 
we can always adjust tq = tq(N) in order to reproduce 
the desired N dependence. 

Let us finally note that apart from the characterstic 
time scale td for translational motion, there are two 
other characterstic time scales in the model, which refer 
to the change of shape and orientation of the ellipsoids. 
These can be identified by the decay of e.g. the correla- 
tion functions (i? 2 ; (t)i? 2 ; (0)) (shape) and (ei(t) • ei(0)) 
(orientation) , where ei is a unit vector pointing in the di- 
rection of the first principal axis (normalized first eigen- 
vector of S a p). An adjustment of these time scales to 
real situations is again possible by tuning Kg properly, 
or by restricting the rotational movements of the ellip- 
soids. For our choice of parameters (K$ = 0.5 and free 
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FIG. 13: Mean square displacement ((r(t) — r(0)) 2 } as a 
function of time for various values of iV in a dense system 
(c = 0.85, e — 1.0). The lines indicate a linear growth at late 
times. Shown for comparison is the mean squared radius of 
gyration (Rq) in the melt for the respective N. 



rotation), td is the largest characteristic time scale. Ac- 
cordingly, for the time scales t>Tr> of interest here, we 
have not to bother about shape-sensitive or orientational 
correlations. 



V. BINARY MIXTURES 

In this section we present the results for a binary mix- 
ture of polymers. We denote the fraction of chains of 
type A with Ja = Ma/M, and the fraction of chains of 
the other type B with fg — Mb/M. The interaction pa- 
rameter is always e = 1.0 and we consider dense systems 
(c = 0.85). 
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FIG. 12: Scaling behavior of the correlation hole. The 
monomer concentration is c = 0.85, the interaction param- 
eter e = 1.0. 



First we determine the coexistence curve. The simula- 
tions are conducted for a semi-grandcanonical ensemble 
0, i.e. we keep the total particle number M fixed and 
allowed the fractions Ja and fg to fluctuate. To this end 
we allow particles to change their identity in the Monte- 
Carlo simulation in addition to the other possible moves. 
The order parameter $ is defined by $=|/a — fg\. For 
T<T C the distribution of the order parameter D(&) has 
a maximum at a finite value $ ma x > 0, while for T>T C , 

*ma X = 0. 

A convenient way to identify the phase transition in the 
simulation is to determine the average ($). For T^lT c , 
D(<I?) becomes sharply peaked and (<&) ~ $ miK . For 
T>T C , ($) is measure of the width of £>(<!>) rather than 
an estimate of $ m ax = 0. However, for T^$>T C the width 
becomes very small and ($) also. Accordingly, the criti- 
cal temperature can be estimated from the region, where 
($) goes from large finite values to zero, see Fig. [14]. 
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We investigated systems of M = 4000, 3000 and 2000 
ellipsoids for chain lengths N = 50, 100 and 200 re- 
spectively, carrying out runs over 4000 MCS. In Fig. |lj 
we show ($) = {\Ja — Jb\) as a function of 1/N5 = 
k-^T /N(cab — e)- The scaling of the coexistence curve 
with N is not perfect since k B T c /S = AN + B [§o|, 
A similar behavior of the phase diagram in dependence 
on N has been found in simulations with explicit chains 
and in experiments [E2| . 



B. Chain dimensions in a blend 

An interesting physical effect in polymer blends is that 
the chains of the minority phase reduce their size to min- 
imize the repulsive interaction with those of the majority 
phase [^0|, [2^| . In Fig. [l5| we present the results for the 
average squared gyration radius of the two different com- 
ponents A and B, divided by the pure melt value (Rq)o 
for <5 = 0. The fraction of B-particles was kept fixed at 
Jb = 0.9. As can be seen from the figure, the particles 
in the minority shrink much more strongly, than those in 
the majority. A detailed discussion of the interplay be- 
tween thermodynamics and chain conformations can be 
found in pil. 



C. Spinodal decomposition 
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FIG. 15: Squared radius of gyration of minority (/a = 0.1) 
and majority (fs = 0.9) type chains as a function of the 
excess interaction parameter S. The simulation was run for 
M = 4000 particles of length N = 50. 



To study the spontaneous phase separation in the GEM 
we consider an equilibrated symmetric blend (/a = 0.5) 
for 5 = and increase the mismatch parameter instanta- 
neously to (5 = 0.2 at time t = 0. The power spectrum of 
the concentration fluctuations is given by the intermedi- 
ate scattering function 



When a polymer mixture is quenched rapidly into the 
mechanically unstable part of the phase diagram, demix- 
ing occurs via spontaneously growing concentration fluc- 
tuations of long wavelengths ]25| . After some initial 
time regime, domains of the two coexisting phases form, 
and with increasing time these domains are expected 
to coarsen according to the Lifshitz-Slyozov growth law 
Lp^^t 1 / 3 , where Ld denotes the average domain size. 



i(M) = (/u(MW(-M)> -4<5(k), 



(39) 
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FIG. 14: The order parameter {\/a — /sj) plotted versus 
l/N5 = k B T/N(i A B - e) for chains of length TV = 50, 100 and 
200. The result indicates, that the coexistence curve scales 
almost as 1/N. 



where pA (k, t) is the Fourier transform of the to- 
tal monomer density g^(y,t) = J^i=i 8a i(y>*) °f A- 
ellipsoids. The power spectrum depends on k= |k| only 
and an explicit expression of I(k,t) in terms of particle 
positions, orientations and sizes is derived in appendix 

The results for I{k,t) are presented in Fig. [l(| At 
early times a peak develops in I(k,t), which is associ- 
ated with a characteristic demixing length. For t < 5000 
MCS the peak position is almost constant in time. At 
later times the domain pattern begins to coarsen and the 
peak of I(k, t) shifts to smaller values of k. To quantify 
this coarsening process we consider the normalized first 
moment 



j^dkkl{k,t) 
f™dkl(k,t) 



(40) 



of the intermediate scattering function, which is a mea- 
sure of the inverse domain size. As shown in Fig. [It] 
the results point to an asymptotic approach towards the 
Lifshitz-Slyozov growth law k\(t) ~ i -1 / 3 . We note that 
in more detailed microscopic models it is hard to reach 
the asymptotic coarsening regime due to the slow diffu- 
sion of the chain molecules (see the discussion in Sec. IV). 



Nevertheless, in some simulations of such models the 
Lifshitz-Slyozov growth law could be seen j26|, p7[. 
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FIG. 16: Intermediate scattering function I(k,t) for various 
times t. In the early time-regime a) the structure develops 
mainly at some fixed characteristic wavelength, while for later 
times b) the peak position of I(k,t) shifts to smaller values 
of k demonstrating coarsening of the domain patterns. The 
function is only evaluated at the physically meaningful integer 
values of k. The data are determined for a system of M = 
4000 particles of size N = 50. 



sequence, no tuning of the interaction parameter e with 
N is necessary to obtain the desired Gauss type behavior 
in dense systems. On the other hand, an important mo- 
tivation of the work of Murat and Kremer for setting 
up the ellipsoid model was to keep a direct link to more 
detailed microscopic models. For such underlying models 
it seems that a tuning of e with N can not be avoided. 
However, one may ask, if it is really necessary to compare 
many particle simulations of the underlying microscopic 
model with many particle simulations of the associated 
ellipsoid model in order to establish the dependence of e 
on N, as it was done in 0. An alternative route, deal- 
ing only with the single chain behavior of the underlying 
microscopic model, would be to use the Gaussian form 
for P(S,N) from eq. (|l5|), where N is the number of 
Kuhn segments of the real microscopic chain with N m { c 
monomers. Then one should include self-interactions (see 
eq. (Q)) with e = e(N) in the associated ellipsoid model in 
such a manner that the true behavior of P(S,iV m i C ) for 
single chains of the microscopic model is best accounted 
for, e.g. by requiring the first moment of Pr(Rq,N) of 
the ellipsoid model with self-interaction to equal the first 
moment of the true Pr{Rq, N m i c ). In this way a later 
tuning of e with ./V in the many particle simulation should 
no longer be necessary. The procedure seems to be partic- 
ularly promising, if effective interactions of longer range 
are included in eq.(Q), which mirror the interactions of 
the monomers in the underlying microscopic model. 

As outlined in the Introduction, our motivation for 
studying the GEM stems largely from our wish to study 
the behavior of thin films of polymer mixtures. In this 
respect we regard the work for the bulk behavior outlined 
here as a suitable basis for further investigation. The im- 
portant influence of substrate interactions on the demix- 
ing kinetics of thin films can be easily taken into account 
by incorporating the coupling of the monomer density to 
an external potential in the free energy functional. Free 



VI. CONCLUDING REMARKS 

In summary we have shown that the GEM pro- 
vides a useful tool to study polymer systems on semi- 
macroscopic time and length scales. The input quanti- 
ties, that are the probability P(S) for an ellipsoid to have 
a shape S and the mean monomer density g(x; S) in the 
principal axis system for a given shape, were investigated 
in detail for Gaussian chains and expressed in terms of 
approximate formulae. Both homogeneous systems and 
binary mixtures have been simulated on the basis of these 
input quantities by means of a discrete Monte Carlo pro- 
cedure. The results from these simulations were shown to 
reproduce the basic scaling relations of polymer physics. 

A valuable feature of the GEM is that the interactions 
between the monomers of the same ellipsoid and different 
ellipsoids can be treated on the same footing by includ- 
ing self-interactions as discussed in Sec. Ill C . As a con- 
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FIG. 17: Time-dependence of the first moment of the interme- 
diate scattering function measuring the inverse domain size in 
double-logarithmic representation. At late stages the asymp- 
totic Lifshitz-Slyozov behavior fci(t) ~ t' 1 ^ 3 is reached. 
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surfaces can be modeled by introducing attractive inter- 
actions between the centers of the ellipsoids. These issues 
will be addressed in the near future. 



Least squares fit of the constants yield u 2 = 0.674 and 
a' 2 = 0.748. 



APPENDIX C: OVERLAP INTEGRALS 
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APPENDIX A: PARAMETERS OF THE 
DISTRIBUTION FUNCTIONS 

The determination of the parameters is based on the 
fact that the moments of the heuristic functions are cal- 
culable analytically. Therefore we define the constants 
a_R and da for Pr(Rq) in such a way, that the first 
two moments (Rq) and (R G ) are exact in the limit 
N — > oo. With this definition and the well known pl| 
results (R%) ~ N/6 and (Rq) ~ 19 iV 2 /540 we are led to 
the two equations 



a.R 



1 K Q (2d R ) K (2d R )K 2 (2dj 
QdnK^dR) ' 



36 Kf(2d f 



19 

540 



(Al) 



At this point we also give a more quantitative test 
of Pr(R g ) by comparison of the higher moments of 
R G , The asymptotic exact result for the third moment 
is (R%)/N 3 = 631/68040 w 0.009274, while we find 
(R%)/N 3 = 0.009280. For the fourth moment the ex- 
act result (R%) /N 4 = 1219/408240 w 0.002986 has to be 
compared with (R%) /N 4 = 0.003000. 
For P a (S a ) the same procedure leads to 



(S a ) K na ^(2d a ) 
d a K„ a -2(2d a ) 
X„ Q _ 3 (2rf a )if na _ 1 (2d a ) 



(S 2 Q 



Kl 



2 {2d a 



(S a ) 



(A2) 
(A3) 



Since no exact analytical result for the moments (S™) 
exists to our knowledge, these quantities are determined 
by Monte-Carlo simulations of Gaussian chains of length 
N = 1000. The results are given in Tab. fv. 



APPENDIX B: IMPROVED FORMULA FOR THE 
MONOMER DENSITY 

For the sake of completeness, we present an improved 
heuristic formula for the density of the x 2 -axis: 



Q2(u 2 ) = 



1 



2tt2ct; 



exp 



exp 



(u 2 - u 2 ) 2 
2 a' 2 

(u 2 + u 2 f 
2 a' 2 



(Bl) 



When inserting eqs.(III B,p8) into cqs.(Q,||) the calcu- 
lation of Fl n tcr reduces to a sum of nine terms of the form 

/+oo 
d 3 yx (1) (y)x (2) (y), (Ci) 
-oo 

where the functions X (y) ? J = 1,2 are given by 



x U) (y) 



7^U) eXP \ 2 



(2^/ 2 nLi^ 

L i) (y) = E^S(^-4' ) )> *$=R$°*- (C2) 



The a a were defined in eq. |lj B| ) and specified in Tab. 0. 



Moreover, = ±K ha R[ 3) u Q + ri j) or = ri j) de- 
pending on which of the three terms in eq. ( p5| ) con- 
tributes to X (y)- We define 




With these definitions the expression (|Cl|) reads 



(2vfa 



+ 00 



dry exp 



T(y) 



(C5) 



Since 



a,/3 



c 



(J) 

a/3 



E 



7?0')-7?(i) 

w 



it follows that 



T(y) = E A a py a yf3 + E B aVa + C 

a,f3 a 



where 



(C6) 

(C7) 

(C8) 
(C9) 





c = EKK'^+^-r)-^ 10 ) 



a,/3 
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The quadratic form (C7) can be rewritten as 

T (y) = Z^iVa - Wot) Da/3 (y ~Wf3)+E, 

which yields 

exp(-5/2) 



(27r) 3 / 2 crVdetP 



4detX> = 4^n^22-433 +^12A3-423 

— Al 3 A22 ~ -4 2 3-4n — ^12^33 



E= C 



1 



16detP 



Bf(A 2 23 ~ 4A 2 iA 33 ) 
Bl(A 2 12 - 4^11^33) 

25x52(2^12^33-^13^23) 
25x53(2^13^22-^12^23) 

2 5 2 5 3 (2^23^ii- ^12^13) 



(Cll) 

(C12) 

(C13) 
(C14) 



Here we defined = Y^ 1 9i ^altfj 3 and g§' 



The center of 

mass of the i'th particle is denoted by r\ . The inter- 
mediate scattering function is then obtained by using 
cq. (D4) in eq. (B9h for all times t (for the A ellipsoids). 



APPENDIX D: FOURIER TRANSFORM OF THE 
MONOMER DENSITY 



For a total monomer density defined by 
1 M 

e ' (y) "M(A7TT)g^ (y) ' 

the Fourier transform is given by 

/oo 
d 3 yg'{y) exp(-ik-y) 
-OG 



(Dl) 



(D2) 



Inserting the q[ (y) from eq. (^) with ( 1IIB 28 ) yields 

P0O - tt^ttX {/></< iky 



2 + ci M 



[x { :Hy) + x^(y) + c lX i l) (y)]}m 



The three different functions x$ (y) > Xh (y) an< ^ X<^ (y) 



are defined in eq. (|C2|), where all different kinds of 



+H 1>a R{ j) u + r { J\ w^ h 



Wa occur, i.e. I»a' a 

— 7£x,e«! Rx Uq + fa and Wa,c — fa\ respectively. The 
integration is straightforward yielding the final result 

M 

P(k) 



1 



1 M 
"a? E 



1 



1 



(2 + ci)cricr 2 cr 3 M ^ y^W v / detO^ 



x exp 



exp 



2 cos Y; k P n il R r u ° + c i 



(D4) 
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